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We report on the search for nearly vertical up-going muon neutrinos from WIMP 
annihilations in the center of the Earth with the AMANDA-BIO detector. The 
whole data sample collected in 1997, 10 9 events, has been analyzed and a final 
sample of 15 up-going events is found in a restricted zenith angular region where 
a signal from WIMP annihilations is expected. A preliminary upper limit at 90% 
confidence level on the annihilation rate of WIMPs in the center of the Earth is 
presented. 



1 Introduction 

There are strong indications for the existence of dark matter in the UniverseE 
The favourite candidate for the dark matter is a Weakly Interacting Massive 
Particle, a WIMP, of which the neutralino that arisesJn supersymmetric ex- 
tensions of the standard model is a natural candidate.H If these WIMPs exist 
they accumulate in the center of the Earth and the Sun, where they can 
annihilate pair-wise producing muon neutrinos that can be searched for with 
neutrino telescopes like AMANDA. Searches of this kind have been performed 
by existing experiments like MACRO 0, Baikal □, Baksan la, Kamiokande and 
Super-Kamiokande B. We here report on searches for muon neutrinos from 
WIMP annihilations in the center of the Earth with the AMANDA-BIO de- 
tector using data from 1997. 



2 The AMANDA-BIO detector 



AMANDA-BIO consists of an array of 302 Optical Modules (OM) arranged in 
10 strings deployed between 1500 and 2000 meters in the South Pole ice cap. 
Muons from charged-current neutrino interactions near the array are detected 
by the Cherenkov light they emit when traversing the ice. The timing of the 
Cherenkov photons reaching the OMs enables us to reconstruct the muon 
track. A more detailed description of the detector is given in U 



3 Signal and background simulations 

Our main backgrounds when searching for a WIMP signal are the atmospheric 
neutrinos and muons. We first discuss the signal and then the background 
simulations. 



3.1 Simulation of WIMP annihilations 

WIMPs annihilate pair- wise to, e.g., leptons, quarks and gauge and Higgs 
bosons. High-energy neutrinos are produced in the decays and/or hadroniza- 
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tion of these annihilation products. Neutrinos produced in quark jets (from 
e.g. bb or Higgs bosons) typically have lower energy than those produced from 
decays of r leptons and gauge bosons. We will refer to the first type of an- 
nihilation channels as 'soft' and to the second as 'hard'. As a typical soft 
channel we choose bb and as a typical hard channel we choose W + W~ above 
the W + W~~ threshold and t + t~ below. The hadronization and/or decay ppf 
the annihilation products were simulated with PYTHIAS For details, see a. 

3.2 Simulation of the atmospheric neutrino flux 

We have simulated the expected atmospheric neutrino flux using the calcula- 
tions of Lipari £3 and the neutrino and anti-neutrino— nucleon cross sections 
from Gandhi et alJiil The actual neutrino-ruicleon interactions have been sim- 
ulated with PYTHIA 13 using the CTEQ3 t3 parameterization of the nucleon 
structure functions. A 3-year equivalent atmospheric neutrino sample with 
energies between lOXJeV and 10 TeV and zenith angles between 90° and 180° 
has been simulated. 113 

3.3 Simulation of the atmospheric muon flux 

The majority of the triggers in AMANDA are induced by muons produced in 
cosmic ray interactions in the atmosphere and reaching the detector depth. 
The simulation of this atmospheric muon flux was performed with an algo- 
rithm described in E-3. Incoming protons were generated between 0° < 8 <85° 
and with energies between 1.3 TeV and 1000 TeV assuming a differential en- 
ergy distribution of E~ 2S7 . Simulating a statistically significant sample of 
atmospheric muon background is an extremely CPU-intensive task due to the 
high rejection factors needed. We have simulated 6.3 x 10 10 primary inter- 
actions, giving 3.9 x 10 7 atmospheric muon triggers, which corresponds to 8 
days of detector lifetime. _ . 

In all these simulations, muon interaction probabilities from tB were used. 

4 Data analysis 

The analysis presented in this paper was performed on data taken with the 10- 
string AMANDA detector between March and November 1997. The original 
data set consists of 10 9 events in a total of 135 days of detector lifetime. We 
apply five different levels of cuts which primarily cut on the zenith angle, 
the number of hits, the track length aryl gome variables that describe the 
space and time topology of the event £3'0 At level 5, we have cut away all 
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Figure 1. Left: Angular distribution of the remaining data (full line) and simulated 
atmospheric neutrino events (dashed line) after filter level 5. The angular range is between 
165° and 180°. The atmospheric neutrino sample has been normalized to the 135 days of 
lifetime of the detector in 1997. Right: Angular distribution of the remaining fraction of 
simulated WIMP signal (mass 250 GeV with hard annihilation channel) at filter level 5. 



the atmospheric muon background and are left with 15 neutrino events in 
data which is consistent with the expected 16.6 events from the atmospheric 
neutrinos. The passing rate for a WIMP with mass 250 GeV and a hard 
annihilation spectrum is 29% at the same level. 

5 Results 

In Fig. [l| we plot the events at level 5 together with the prediction for the 
atmospheric neutrinos .113 We also plot the expected angular distribution from 
a WIMP. As can be clearly seen, we do not see any statistically significant 
excess of nearly vertical neutrino-induced muons. Hence, we will use these 
remaining events to derive an upper limit on the WIMP signal. 

5.1 Limits on the WIMP signal 

As seen in Fig. |l], the WIMP signal is very peaked towards vertical muons, 
suggesting that we should cut further in the zenith angle. If we cut at 6 = 
165°, 170°, 172.5° or 175° we can derive 90% C.L. limits on the number of 
signal events, 7.0, 5.6, 4.5 and 5.2 events respectively. Although the systematic 
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Table 1. The 90% confidence level on the neutrino to muon conversion rate (r„ M ) (for 
muons-iabove 10 GeV), the WIMP annihilation rate (r^) and the muon flux (above 1 
GeV)H. 



uncertainties for different angular cuts currently are under investigation, we 
choose here to use the angular window which gives the best limits. From these 
limits on the number of signal events we can calculate limits on the neutrino- 
to-muon conversion rate, r M „, which is the number of neutrinos that interact 
and produce a muon per volume element per second. In Table [l] we show the 
limits on the neutrino-to-muon conversion rate where the muon energy (at the 
neutrino-nucleon vertex) is above 10 GeV. We also show the corresponding 
limits on the WIMP annihilation rate in the center of the Earth, Ta and to 
compare with previous experiments, we also show the calculated limits on the 
neutrino-induced muon flux above 1 GeV. 

In Fig. ||, the limits on the annihilation rate are shown and in Fig. || we 
show the limits on the neutrino-induced muon flux. In that figure we also show 
predictions from the MSSM □ when the WIMP is a neutralino and we also 
compare with other experiments. We see that, especially at higher masses, 
AMANDA is already with this limited data set competitive with previous 
experiments. 

5.2 Systematic uncertainties 

The limits are subject to experimental and theoretical systematic uncertainties 
which affect the estimated background and the effective volume and propa- 
gate to our derived limits. The precise effect is currently under detailed in- 
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Figure 2. 90% confidence level upper limits 
on the WIMP annihilation rate in the center 
of the Earth, as a function of the WIMP 
mass and for the two extreme annihilation 
channels considered in the analysis. 



Figure 3. 90% confidence level upper lim- 
its on the neutrino-induced muon flux for 
a hard annihilation spectrum. The limit 
on the muon flux has been calculated from 
the limits on the annihilation rate with a 
muon energy threshold of 1 GeV to be able 
to compare with other experiments. Also 
shown are predictions from calculations in 
the MSSM.EJ 



vestigation but preliminary estimates indicate that the limits may have to be 
increased by about a factor of 2-3. Ell 



6 Conclusions 



We have performed a search for an excess of vertically upgoing muons with the 
AMANDA-B10 neutrino detector as a signature for WIMP annihilations in 
the center of the Earth. Limits on the signal from WIMPs have been derived 
from the non-observation of an effect. The limits presented here are obtained 
with 135 days of effective detector lifetime and are already comparable with 
limits obtained by Baikal, Baksan, Super-Kamiokande and MACRO, with 
much longer accumulated lifetimes. The systematic uncertainties are under 
investigation and preliminary results suggest that the limits may increase by 
about a factor of 2-3 when they are included. 
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